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The presence of weak intergalactic magnetic fields can be studied by their effect on electro- 
magnetic cascades induced by multi-TeV 7-rays in the cosmic radiation background. Small deflec- 
tions of secondary electrons and positrons as the cascade develops extend the apparent size of the 
emission region of distant TeV 7-ray sources. These 7-ray halos can be resolvable in imaging at- 
mospheric Cherenkov telescopes and serve as a measure of the intergalactic magnetic field strength 
and coherence length. We present a method of calculating the 7-ray halo for isotropically emitting 
sources by treating magnetic deflections in the cascade as a diffusion process. With this ansatz the 
moments of the halo follow from a set of simple diffusion-cascade equations. The reconstruction 
of the angular distribution is then equivalent to a classical moment problem. We present a simple 
solution using Pade approximations of the moment's generating function. 

PACS numbers: 95.85.Pw, 98.62.En, 98.70.Rz, 98.80.Es 



I. INTRODUCTION 

The presence of large-scale magnetic fields in cosmic environments can be probed by various astronomical techniques. 
Synchrotron radiation of relativistic electrons can be detected by its characteristic linear polarization and spectrum. 
Faraday rotation of linearly polarized emission tests the birefringent properties of a dilute magnetized plasma filling 
intergalactic space. Zeeman splitting of an atom's energy levels can be observed by the corresponding shift of spectral 
lines from astrophysical masers. With these standard methods it has been possible to identify micro-Gauss magnetic 
fields coherent over galactic scales in many galaxies and galaxy clusters [2 [2] . 

The origin of these large-scale magnetic fields is unclear. It is assumed that galactic magnetic fields can be main- 
tained and amplified via a dynamo mechanism, where the kinetic energy of a turbulent interstellar plasma is converted 
into magnetic energy [3j. However, this requires initial seed fields of unknown origin, possibly pre-galactic or primor- 
dial [H [5] . The strength and correlation length of primordial intergalactic magnetic fields (IGMFs) can be limited by 
their effect on various stages in cosmic history. The strongest bounds on the strength of primordial IGMFs arise from 
the study of temperature anisotropics in the cosmic microwave background (CMB) [S]. The limits are of the order 
of nano-Gauss for a correlation lengths larger than a few Mpc. Simulations of large-scale structure formation favor 
long-range IGMF with a strength of the order of pico-Gauss [7]. 

It has been suggested that weak IGMFs of the order of fempto-Gauss can be probed by their efi^ect on electro- 
magnetic cascades initiated by the emission of distant multi-TeV 7-ray sources [HI E] ■ High-energy 7-rays produce 
pairs of electrons/positrons in the cosmic infrared/optical background (GIB) with an interaction length of the order 
of 100 Mpc. The secondary leptons lose their energy via inverse-Compton scattering off the background photons and 
produce secondary 7-rays at somewhat lower energies. If these photons are still above the pair-production threshold 
the cycle repeats. In this way the electro-magnetic energy of the cascade is continuously shifted into the GeV-TeV 
energy region. In the presence of magnetic fields secondary leptons are deflected off the line-of-sight and secondary 
7-rays inherit this deflection. This will attenuate the flux originally emitted towards the observer. However, 7-rays 
initially emitted away from the observer can be scattered back into the line-of-sight and partially compensate for this 
loss. 

There are various ways to infer the strength Bq and correlation length As of the IGMF from this effect. For 
small deflections and isotropically emitting sources (or sufficiently large jet opening-angles) the net effect will be an 
extended emission region of secondary 7-rays [51 [TDHT5] . For burst-like 7-ray sources this can also cause an observable 
time-delay between the primary burst and secondary 7-rays [9l [14] . In the case of a hard TeV 7-ray emission the 
secondary component can dominate over the attenuated primary 7-rays. Non-observation of the point-source in the 
GeV-TeV band can then imply a lower limit on the magnetic field depending on the instrument's resolution [5]. These 
methods have been applied to various TeV 7-ray sources [T5H^ and indicate the presence of an IGMF. The inferred 
lower limits on its strength range from 10~^^ G to 10""'^^ G, depending on many systematic uncertainties like the 
primary emission spectrum, the GIB and the coherence length of magnetic fields. 

Besides the systematic uncertainties of these methods, there are also some technical challenges in calculating the 
energy and angular spectrum of the 7-ray halos. A straightforward Monte-Garlo calculation of the electro-magnetic 
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cascade can become numerically expensive; since energy is conserved in the cascade the number of 7-rays, electrons 
and positrons in the cascade increases by one order of magnitude for every decade of the energy shift. At every step 
of the cascade each particle will have accumulated a deflection angle with respect to the line-of-sight which depends 
on its history in the cascade. In order to accumulate a satisfactory resolution in energy and angular extend of the 
halos it is necessary to sample over many cascades. 

In the absence of deflections by magnetic fields the electro-magnetic cascade can be calculated efficiently by ana- 
lytical methods using cascade equations and the method of matrix doubling [55]. We will show in this paper that 
there is a straightforward extension of this method to the case of isotropic emitters and small deflections in magnetic 
fields. The key observation is that the deflection of electrons and positrons in the cascade in combination with 
inelastic losses to photons can be treated as a diffusion process in 0-space where the diffusion coefficient depend on 
the particle's Larmor radius and the inverse-Compton energy loss length. We derive diffusion-cascade equations that 
describe the evolution of the moments of the 6'-distribution and give a simple method how these moments can be used 
to reconstruct the distribution. 

We will begin in secti on [Tl| by a discussion of electro-magnetic cascades from 7-ray point sources in the presence of 
weak IGMFs. In section |III| we will motivate the extension of the Boltzmann equations by a diffusion term in 0-space 
and give an extended set of cascade equations for the moments of the ^-distribution. We discuss in section |IV] how 
the full ^-distribution can be reconstructed efficiently from a finite number of moments via explicit inverse Laplace 
transformations of Fade approximations of the moment's generating function. We will test our method in section [V] 
by two examples and compare our results to previous studies. Finally, we conclude in section [Vl} 

We work throughout in natural Heaviside-Lorentz units with /i = c = eg = /io = 1, a = e^/(47r) ~ 1/137 and 
1 G ~ 1.95 X 10-2eV^. 



II. ELECTRO-MAGNETIC CASCADES 



The driving processes of the electro-magnetic cascade in the cosmic radiation background (CRB) are inverse Comp- 
ton scattering (ICS) with CMB photons, e^-|-7bgr 6^+7, and pair production (PF) with CMB and GIB radiation, 
7 -I- 7bgr — + e~ [231 El]- In particular, the interaction length of multi-TeV 7-rays depend on the GIB back- 
ground at low redshift and is of the order of a few 100 Mpc. We show the relevant interaction lengths and energy 
loss lengths in the left panel of Fig. [T] High energetic electrons and positrons may also lose energy via synchrotron 
radiation in the intergalactic magnetic field, but this contribution is in general negligible for the small magnetic field 
strength considered here. Further processes contributing to the electro-magnetic cascade are double pair production, 
7 + 7bgr — >■ e+ -|- e~ -I- 6+ -I- e~ , and triple pair production, -I- 7bgr e"^ + + e~ [221 US]. These contributions 
can be neglected for cascades initiated by multi-TeV 7-rays considered here. Also, interactions on the cosmic radio 
background are negligible in this case. 

For the calculation of the flux from a 7-ray point-source it is convenient to start from the evolution of a comoving 
number density Yq — nal{\ + z)^ (GeV~^ cm"-'^) in a spatially homogeneous and isotropic universe. The Boltzmann 
equations of electrons/positrons (Ye) and 7-rays (Ky) is given by 

Y^{E) = dE{HEY^)-T^Y^{E)+ J dE'-fp^{E' , E)YpiE') + C^{E) , (1) 

together with the Friedman-Lemaitre equations describing the cosmic expansion rate H{z) as a function of the 
redshift z, H'^(z) = Hq [ftm{^ + z)'^ + ^a], normalized to its present value of Hq ^ 70 kms^^ Mpc~^. We consider 
the usual "concordance model" dominated by a cosmological constant with fl^ ^ 0.7 and a (cold) matter component, 
ilm ~ 0.3 [2^. The time-dependence of the redshift is given by dz — — dt(l -I- z)H. The first term in the r.h.s. of 
Eq. ([1]) accounts for the continuous energy loss due to the adiabatic expansion of the Universe. The second and third 
terms describe the interactions with background photon fields involving particle losses (a — >■ anything) and particle 
generation a — >■ /3. The angular-averaged (differential) interaction rate, Ta (lafs) is defined as 

1 

Ta{z,Ec,) = ^ J dcos6' J de (1-/3 cos 6')n^(z,e)cr^°^*, (2) 
-1 

^^piz,E^,Ep)^Ta,iz,Ea,) ^^{Ea,,Ep), (3) 
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where n^{z, e) is the energy distribution of background photons at redshift z and dN^p /AEp is the angular-averaged 
distribution of particles (3 after interaction of a particle a. Besides the contribution of the CMB we use the CIB from 
Rcf. [37]. Due to the cosmic evolution of the radiation background density the interaction rates ^ and ([S]) scale 
with redshift. The CMB evolution follows an adiabatic expansion, n^{z, e) = (1 + z)^ n^(0, e/(l + z)), and we assume 
the same evolution of the CIB for simplicity. We refer to Ref. [5S1 for a list of the redshift scaling relations of the 
interaction rates in Eqs. ([T]). The last term in Eq. ([I]), accounts for the emission rate of particles a per comoving 
volume. 

In the limit of small deflections of particles via magnetic fields, the flux from a 7-ray point source at redshift distance 
z'* with emission rate is equivalent to an angular-averaged flux from a sphere at redshift z* . Hence, the solution 
of y at t = is equivalent to the point source flux J (GeV^^ cm^^ s~^) by replacing the emission rate density C in 
by 

where the comoving distance of the source (in a flat universe) is given by dc{z) = J^dz'/H{z'). Note, that we can 
also use the ansatz Q for a cosmic ray (CR) point source located at redshift z*, where the electro-magnetic emission 
is in the form of cosmogenic 7-rays, electrons and positrons produced during CR propagation |29j . 



III. ANGULAR DIFFUSION IN INTERGALACTIC MAGNETIC FIELDS 



The 7-ray cascade can only contribute to a GeV-TeV point-source flux if the deflections of secondaries off the 
line-of-sight are sufficiently low. The scattering angle of secondaries is only of the order of e/we and can be neglected 
for the optical/infra-red background photon energies e. However, electrons and positrons can be deflected in the 
IGMF. We can estimate the extend of the cascaded 7-ray emission by simple geometric arguments following [llj . 
Deflection of electrons and positrons will be small if the energy loss length Ag of electrons/positrons via inverse 
Compton scattering (ICS) is much smaller than the Larmor radius given as Rl = E/eB ~ l.l(£'Tcv/^fG)Mpc. Here 
and in the following we use the abbreviations E = i^xcvTeV, etc. For center of mass energies much lower than the 
electron mass, corresponding to energies below PeV in the CMB frame, electrons and positrons interact quickly on kpc 
scales but with low inelasticity proportional to their energy, Ag — 0.4 Mpc/i^xcV- The typical size of the point-spread 
function (PSF) of imaging atmospheric Cherenkov telescopes (lACTs) is of the order of Spgp — 0.1°. Hence, magnetic 
deflections become important if 6'psf < Xe/RL or E<U TeVy/B^ 

PSF,0.1°- 

In the following we are going to study these magnetic deflections more quantitatively. For simplicity, we will start 
with a regular IGMF that fills the space between the source and the observer and has the component B± perpendicular 
to the line-of-sight. We also assume that the source is emitting 7-rays isotropically. ^ Due to charge conservation in 
the cascade electrons and positrons will be produced in equal rates and will be deflected in opposite directions. For 
small scattering and isotropic emission we can assume that leptons that are lost by deflections out of the line-of-sight 
are replenished by the corresponding leptons deflected into the line-of-sight. Effectively, we can hence assume that the 
total number of electrons/positrons within the line-of-sight remains constant by these deflections while the scattering 
angle 9 is broadened by the magnetic field. The width of this 6'-distribution, ye{E,6), is determined by the energy 
loss length via ICS. Secondary 7-rays will inherit the 6'-distribution of the parent leptons and will appear as extended 
halos. 

The energy loss length via ICS with CMB photons is much smaller than the typical distance of TeV 7-ray sources 
or the interaction length of PP in the CIB. This indicates that we can treat magnetic deflections in the cascade as 
a diffusive process of the angle 9. The mean free path of the electrons/positrons corresponds to the energy loss rate 
in ICS and the diffusion velocity is the inverse Larmor radius. Hence, the diffusion coefficient D is of the order of 
AjQg/i?!^. A more rigorous derivation (see Appendix IaI) shows that the evolution of the 6'-distributions of leptons and 



We can relax this condition by requiring that the 7-ray emission is into a jet with sufficiently large jet opening-angle. 
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7-rays, y^E, 9) and y-y{E, 9) respectively, can be described by the coupled set of differential equations, 

~ dE{HEy^) - T^y^ + J2 J <^e'^^^{e' ,E)y^{E') + c*^5{e) , (5) 

Q=e,7 E 

3>e ~ dE{HEy,) - T,y, + J2 [ d^'7oe(£^', E)y^iE') + Cl5{e) + / dE'V,,^{E', E)dly,{E') . (6) 
The diffusion matrix of electrons/positrons in a regular magnetic field is given by 

= ET.cm E'^ 

where {x){E) is the inelasticity of ICS with interaction rate rics(£^)- For cosmological sources the redshift scaling of 
the diffusion matrix (7| can also become important. For primordial magnetic fields scaling as B^{z) = (1 + z)^i?j^(0) 
and ICS with CMB photons the redshift dependence is given by the simple relation I?i.eg(2:, E' , E) — (l + z)'^I?reg(0, (1 + 
z)E',{l + z)E). 

This formalism has the advantage that we can calculate the moments of the ^-distribution by an extended set of 
cascade equations as we will see in the following. Firstly, we introduce the quantities 

oo 

^^h^W)\ I ^^^'"^^/^ (regular). (8) 

— CO 

At leading order we have F^^-* = F as the solution of Eq. ([I]) and for n > 1 the quantities F^"^ correspond to the 
scaled moments of the ^-distribution.^ It is easy to see that the quantities F^"' [n > 0) follow the coupled set of 
differential equations, 

ri")(£;) = aB(Hi?ri"')-rari")(£^)+ E j '^e'^p^{e',e)y^''\e') + 6,„ j dE'v{E',E)Yi^-^\E'), (9) 

0=en E E 

in addition to Eqs. Q. Note, that electro-magnetic interactions of photons and leptons that drive the cascade happen 
on time-scales much shorter than the rate of adiabatic losses in the Universe. We can hence treal-the interaction rates 
as constant over small time- intervals and neglect the energy loss terms dE{HEY^^) in Eqs. (Ilj) and (oj). We show 
in Appendix [B] that this system of equations can then be solved efficiently by a generalization of the conventional 
cascade equations. 

We next consider a randomly oriented IGMF field with a coherence length As much smaller then the distance d of 
the source. In this case we have to replace Eq. ([5| by the evolution of radial diffusion on a sphere of the form'^ 

y,^dE{HEy,)-T,y,+ ^ J dE'-f^e{E',E)y^{E') + C;5{9) + J dE'V,^d{E',E)e-'de[9deye{E')] , (10) 
with diffusion coefficient (see Appendix [A| 

l min(l,ABrics(i^)) e^Bg 
^"'^(^ '^^ - 3 ET.csiE) i?'2(x)(i?')- ^ ^ 

Here, a factor 1/3 accounts for the random orientation of the magnetic field w.r.t. the line-of-sight. Analogously to 
the diffusion in a regular magnetic field we can define moments of the diffusion in random IGMFs by 



Y 



in) 



2tt 



h (2"n!)2 



CX3 

Jd9 99^''ye/j (random), (12) 



^ To be more precise, 9 is an element of the covering space R of the circle S^. The distribution along the circle is then obtained by 
^gi {E, 9) = X]nGZ ^rC^i ^ + 27rn). However, we are only interested in small scattering angles 9 <^ 1° and hence y^i {E, 9) ~ yuiE, 9). 
^ We consider only small deflections and can hence approximate the sphere as two-dimensional flat space. 
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FIG. 1: Left panel: The interaction length (sohd lines) and energy loss length (dashed lines) from various contributions of 
the CRB. We show the rates separately for the CMB and GIB. Also shown is the inverse Hubble scale (dotted line). Right 
panel: The spectra of 7-rays from a source at 120 Mpc with injection spectrum ~ exp(— _E/300 TeV) (gray line) 
following Ref. |12| . We show the contribution of surviving primary 7-rays (dashed line) and secondary cascaded 7-rays (solid 
line) separately. 



which follow the same differential equations ([9| with diffusion matrix I?rnd- 

So far we have only considered the diffuse scattering of the photons along their initial trajectory. How does this 
translate into the observed morphology of the 7-ray signal? Deflections of electrons close to the source at distance d, 
e.g. by the magnetic field of the source itself, will have a weaker impact on the observed angular distribution than 
deflections close to the observer. If the cascade experiences a deflection /S.9 at a distance r from the observer we can 
approximate the corresponding angular displacement A0' in the observer's frame via A0'/A0 ~ (d — r)/d. We can 
account for this scaling in the cascade equation by introducing the corresponding scaling in the convection velocity 
R~j^^ or, equivalently, by a scaling of the diffusion matrix of the form V ~ ((d — r)/d)'^'D. In practice, this requires 
that we repeat the calculation of transfer matrices after sufficiently small propagation distances, for which we then 
also account for the variation of (differential) interaction rates F (7) with redshift and adiabatic energy loss. With 
this simple modification the moments F'-"-' reflect the angular distribution of 7-ray halos, as long as scattering in the 
magnetic field is small and the source is emitting isotropically. 

As an example, we study in the following an isotropic 7-ray point-source at a distance of about 120 Mpc - as Mrk 421 
- with a 7-ray injection spectrum of the form ^ exp(— i?/300 TeV). This particular example has been studied 
in Ref. |12) and hence our results are directly comparable. In the left panel of Fig. [l] we show the source spectrum, 
i.e. the spectrum that would be visible without the CRB (thin gray line) together with the electron/positron and 
7-ray spectrum after propagation. The total 7-ray spectrum (dotted green line) can be decomposed into a "primary" 
component (dashed green line) of surviving 7-rays and a "cascaded" component (solid green line) from 7-rays of the 
cascade. The 7-ray flux is strongly suppressed beyond 10 TeV due to the PP with the CIB and secondary 7-rays from 
ICS with the CMB peak between 0.1-1 TeV. 

We will assume in the following that the cascade develops in a weak IGMF with strength ~ 10"^^ G and a 
coherence length — 1 Mpc extends. For the reconstruction of the 7-ray halo it is convenient to first subtract the 
moments of the surviving primary 7-rays that don't take part in the cascade, 

7, halo 7, total 7 primary ' V 

In our example we assume a point-source with sufhciently small angular extend, corresponding to the case X^^'^pj-imary ~ 
Yry^ primary and vanishing higher moments. In general, the higher moments of the primary source with an angular extend 
29a can be approximated by 

primary P'^""'*''^ n' (tZ -|- 1) '4" ' 
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The size of the first non-trivial moment Xy^^haio/^7°haio aheady serves as a first indicator for the size of the 7-ray 
halo. If this is much larger than the PSF of an lACT the flux of secondary 7-rays will be strongly isotropized and 
can only be constrained by the diffuse 7-ray background (see e.g. [SO])- We will show in the following that we can use 
the spectrum of moments to reconstruct the 7-ray halo for small deflection angles. This will also give an indication 
at which energies the contribution of secondary 7-rays contribute to the point-source spectrum. 



IV. RECONSTRUCTION OF THE ANGULAR DISTRIBUTION 



The moments of the 7-ray halo serve as a measure for its angular distribution. How we can reconstruct the angular 
distribution from a limited number of moments? As a first step it is convenient to define a distribution f{E,x) by 
the transformation 



3^7,halo(£^,^) = i^7,halo(£^) / dx 





1 



(27rx) 



-e'^/2x 



f{E,x)., 



(15) 



for regular {d — 1) or random (d — 2) magnetic fields. This transformation is motivated by the observation that the ker- 
nel Gd{x, 9) — [. . .] corresponds to a Green's function of the d-dimensional diffusion equation, {dx — ^^ dg.)Gd{x, 9) ~ 

and Gd{0,9) = JJ^S{9i) with J2i^i — We can then identify the quantities y(")/F(") as (scaled) moments of the 
distribution f{E,x) for both, regular and random fields: 



00 

lin{E) = J dxx'' f{E,x)^r 



(16) 



We hence arrive at a classical (Stieltjes) moment problem |31l of finding the distribution / from its moments 
From the differential equations ^ and the definition (16 1 it is easy to see that we can find a constant C 
such that finiE) < Cn\[2dmaxE'>Ei'DXE' , E))y^ , where d the distance to the source. This is a sufficient condition 
for a determinate moment problem (31j . i.e. there exists a unique solution / satisfying Eq. (16). Note, that the 
reconstruction of / from the complete set of moments /j,„ is trivial. For instance, we can express / by an infinite 
sum of Laguerre polynomials which are orthogonal on [0, 00) under the measure exp(— x). However, this method does 
not prove convenient if there are only a finite number of at our disposal. The truncation of the expansion after 
the first -f- 1 basis function leads typically to rapidly oscillating solutions. Alternatively, we can reconstruct the 
distribution by a sequence of approximations /, which are maxima of an entropy functional j32j . where the condition 
( 16 1 are introduced via Lagrange multipliers. This problem can then be reduced to a minimization problem of an 



A^-dimensional effective potential. 

In our case we choose a different approach, which is suitable for the particular form of the distribution. First, we 
introduce the Laplace transform of the potential / as 



f{E,s)^C{fiE,x)} 



00 

/OO 
dxe-^^fiE,x)^Y. 

n fe=0 



i-sl 

kl 



(17) 



The Laplace transform / corresponds to a generating function of the moments, (— l)"9"/(i?, s)|s=o = Mri(-B), and the 
solution to the moment problem corresponds to the inverse Laplace transform f{E,x) — £^^{f{E,s)}. However, in 
practice we have only a finite number of moments A^ + 1 and the truncation of the alternating scries ( |17[ ) does not 
converge for large s. 

We can find an approximate solution by replacing the truncated series by a Fade approximation - a method which 
is well-known to chemistry, engineering or nuclear physics [33] . We are approximating / by a rational function 
f[M,M+i]{s) — P{s)/Q{s) where P and Q are polynomials of degree M and M + 1, respectively. The coefficents of P 
and Q are determined by matching the first 2M ~|- 1 terms of the Taylor expansion of f[M,M+i] to the truncated series. 
Clearly, for A^+ 1 calculated moments we can only consider M < N/2 for the approximation. Since deg(Q) > deg(P) 
the Fade approximation is finite as s — ^ cxd, in contrast to the truncated series it approximates. If we write the 
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FIG. 2: Left panel: The first 15 non-trivial moments fj,„ at three different 7-ray energies for the example shown in the right 
panel of Fig. nland assuming a IGMF with strength Bo — 10-^= G and A_b = 1 Mpc. The dashed lines show the momenta 
reconstructed by the Fade approximations /[4,5] that are fixed by the first 6 calculated moments. The approximation reproduces 
the higher moments well. Right panel: The 7-ray halos reconstructed from the moments shown in left panel. We also indicate 
the typical size of the PSF for lACTs. 



denominator via its roots Si with multiplicity rrii 
rational function /[m,m+i] has the simple form 



Q{s) — nr=i('5 ~ the inverse Laplace transform of the 



f{E,x) ~ £ H/[Af,M+i] 



} = EE 



,iE) 



=1 j=i 



ij 1)! 



(18) 



where the coefficients follow from an expansion into partial fraction. Note, however, that for a general Pade 



approximation it is not guaranteed that all 5R(si) < and hence the approximation (18) can be unstable even if the 



exact solution (17 1 is stable itself. However, by lowering the degree of approximation M it is in general possible to 
obtain a stable Pade approximation that fulfills the necessary criteria. This can be done by trial and error - as we do 
here for simplicity - or by an algorithmic procedure [33]. We will show in the following section, that this procedure 
is stable and reproduces the moments of the distribution well. Finally, the distribution N(9) can be obtained from 
Eqs. and (flSl). 



We illustrate this procedure for the cascade spectrum shown in the right panel Fig.[T] In the left panel of Fig. [2] we 
show the first 15 non-trivial moments of the distribution / for 7-ray halos at 10^'^, 10'^ and lO'^'^ GeV. The dashed 
line shows the moments calculated via the Pade approximation /[4^5] . Note, that this approximation is determined by 
the first eight non-trivial moments, but also reproduces all the higher moments of our calculation satisfactorily. Using 
Eqs. ( jTs] ) and ( 18 1 we can derive the angular distribution of the halos which are shown in the right panel of Fig. [2] 
For illustration we normalize the distribution as N{9,E) = 3^7,haio(^j £')/^, total (£')■ Note, that not all of this 7-ray 
halo will be resolvable in lACTs. We are indicating in the plot the typical size of the PSF of 0.1°. We will discuss in 
the following section the size of these 7-ray halos in more detail. 



V. SIZE OF THE EXTENDED HALOS 



The size of the extended halo serves as a measure of the IGMF. Typically, the low-0 form of the halo derived from 
the approximation (18 1 depend on a few roots with large real component |3fi(si)|. In this case, the ^-distribution is 
in the form of a modified Bessel function for a random IGMF with Xb d. The sub- halos have the form 



1 



1 



(19) 
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FIG. 3; Left panel: The size of the extended 7-ray halo defined by Eqs. ( 19 1 and ( 20 \ for a source a.t z — 0.031 with spectrum 
Q-yiE) ~ £"^exp(-£;/300TeV). Right panel: A model for the 7-ray spectrum of the blazar source 1ES0229+200 located 
at 2; = 0.14. The blue data points show the H.E.S.S. observation [35] and the red lines correspond to the upper flux limits 
from Fermi-LAT (taken from Ref. [18]). We assume a source spectrum of the form oc i5~^''^0(2OTeV — E). The solid 
green line shows the spectrum of secondary 7-rays without deflections in the IGMF. The dotted green lines indicate the part 
of the cascaded 7-ray spectrum within 0.1° around the source for an IGMF with coherence length Xb ~ 1 Mpc and strength 
Bo = lO"^** G, 10"^^ G and 10"^* G, respectively. 



where in the last step we took the asymptotic form of Kq at large 9 and introduced the characteristic size of the 
sub-halo, 9i = (2|5R(s,)|)-i/2. Hence, there is a simple relation between the measurable size of the halo and the simple 
poles of the Pade approximation. 

In the cases shown in the right panel of Fig.[2]the leading order halo is below the typical instrument's resolution of 
^PSF = 0.1°. Instead, the next-to- leading-order root will determine the size of the halo. In general, we hence define 
the size of the leading (observable) halo as 

0,^t = mm{{9,}\e, > 9psf) • (20) 

For our test spectrum we show the parameter ^cut in the left panel of Fig. [3] for various magnetic field strengths and 
7-ray energies between 100 GeV to a few TeV. As before we consider a coherence length of Ab = 1 Mpc. The size 
of the halo in this energy range follows approximately ^cut oc as the fit shows. This agrees with the findings of 
Ref. P2] (Fig. 7) derived from a Monte-Carlo study. 

Another interesting situation occurs if the cascaded 7-ray spectrum dominates over the primary 7-ray emission. 
This can happen for injection spectra that are considerably harder than E~'^. In this case the detection sensitivity of 
the cascaded GeV-TeV spectrum depends on the size of the halo and the resolution of the telescope. As an example 
we consider the emission of the blazar source 1ES0229-I-200 located at redshift z = 0.14, which has been detected by 
its TeV 7-ray emission by H.E.S.S. [35]. The spectrum is shown in the right panel of Fig.|3]as the blue data. Following 
Ref. [TS] we model the 7-ray emission spectrum as Qj oc i?^^/'^9(20TeV — E) (thin gray line). The surviving primary 
7-rays are shown as a dashed green line and secondary cascaded 7-rays by a solid line. 

It is easy to understand the shape of the various spectra. Primary 7-rays close to i^max interact with the GIB to 
produce electron/positron pairs. This is a slow process happening on typical scales of the order of a few 100 Mpc (see 
left panel of Fig.Q. The lept ons quickly lose energy via ICS with CMB photons at a rate bics — E/X^; their spectrum 
in quasi-equilibrium (dtYe ~ 0) follows the differential equation 9_B(6ics^e) — PppKy. Thus, the Comptonized electron 
spectrum for E ^ i?max has the form ~ E~'^. The typical photon energy from ICS of a background photon with 
energy e is given by i?^ ~ e{Ee/me)^ . The resulting photon spectrum at _E ^ ^-max follows from energy conservation 

in ICS, dtY^ ~ {AEe/AE^){bics/E^)Y^ - {E^/ E^f / {2\e)Ye - E~^''^ . The plateau of the full cascaded spectrum 
shown in the right panel of Fig. |3] is slightly softer than this since a part of the inverse-Gompton spectrum is still 
above pair-production threshold and enters a second cascade cycle. 

We also show the expected contribution of secondary 7-rays confined within the PSF of a typical lACT with 
^PSF = 0.1° assuming an IGMF with coherence length As = 1 Mpc and strength Bq — 10^^^, Bq — 10^^^ and 
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Bq = 10 G, respectively. The deflection of an electron of the Comptonized spectrum is approximately 9e ~ Xe/Rh/^ 
following from dt{9eYc,) ~ Y^/Rl and Ye ~ . This is consistent with the results of our diffusion ansatz since the 
first moment of the electron/positron distribution follows Fi^V^e — ('^e/^L/4)^/2. If the typical deflection 9^ 
exceeds ^psf we expect to see a reduction in the point-source flux by a geometric factor (^psp/^e)^ oc E'^^ ^ -^7 ^• 
For ICS in the CMB the transition is expected to occur close to the energy 



This agrees well with the reduced cascade flux ("0 < 0.1°"; dotted lines) shown in the plot. 

Before we conclude we would like to emphasize a subtlety concerning the contribution of the CIB in ICS. As can be 
seen from the summary of interaction/loss lengths in the left panel of Fig. [l] the contribution of the CIB to the total 
energy loss of ICS is negligible. The 7-ray spectrum is hence almost independent of this contribution, but this is 

not the case for the higher moments Ylj"'\ To see this, let us consider a fully Comptonized electron/positron spectrum 
Yg ~ E~'^. Following our previous arguments we have dt{9jYj) ~ {dEe/dEj){bics/ Ej)9eYe ~ [E^/ E^Y / (9,RL)Ye for 
the Comptonized electron spectrum. The growth of the deflection is hence proportional to \/e and optical photons are 
expected to contribute much stronger than CMB photons. However, the fraction of photons that contribute with this 
large deflection is negligible. Inverse-Compton scattering by the CIB will form a shallow plateau of 7-rays that are 
negligible for the calculation of the moments of the central halo from the CMB contribution. We can hence neglect 
this contribution in the calculation of moments which improves the quality of the halo reconstruction at low 9. 



We have discussed a novel technique of calculating extended halos of TeV 7-ray sources in the presence of inter- 
galactic magnetic flelds. The method builds on standard cascade equations that account for all particle interactions 
with the background radiation and treats the effect of secondary electron/positron deflections in intervening magnetic 
fields by a diffusion ansatz. The moments of the angular distribution can be calculated efficiently by an extended set 
of cascade equations. The first moments of the distribution already serve as a good estimator of the halo size. We have 
shown how the full distribution can be reconstructed from further moments via an inverse Laplace transformation of 
the moment's generating function using Fade approximations. 

Our method applies to situations where the emission of 7-rays is isotropic or within a sufficiently large jet opening- 
angle. The 7-ray halo is expected to show further structure in the more general case. For instance, 7-ray emission 
into narrow jets are expected to produce additional breaks in the halo profile TP and non-spherical geometries in the 
case of an off-axes emission |36 . For the illustration of the method we have considered a steady 7-ray emission. In 
the case of pulsed or short-lived 7-ray sources there will be a time-delay between primary and secondary 7-rays due 
to the increased path length of the leptons. This can also serve as a measure for the intergalactic magnetic field. 



This work is supported by US National Science Foundation Grant No PIIY-0969739 and by the Research Foundation 
of SUNY at Stony Brook. 



We assume in the following that the magnetic field is perpendicular to the line-of-sight of the source. In this setup 
electrons and positrons will be deflected by an angle in a plane normal to the magnetic field. For isotropic emission 
and small deflections the leptons deflected off the line-of-sight are replenished by leptons initially streaming away from 
the observer. The net effect is a broadening of the 6'-distribution {y±) due to a convection term with opposite sign 




(21) 



VI. CONCLUSION 
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Appendix A: Derivation of Eqs. ([s]) and ( 10 1 
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for electrons (— ) and positrons (+). The transport equations take the form 



3>± {E, t,9)^± --l—dgy± {E, t,9)-T, {E)y± {E, t, 9) 
^l[E) 



+ J dE'(^^j^,{E',E)y^{E',t,9)+j,,{E',E)y±{E',t,9)^ +^C:{E,t,9 

E 

The evolution equation of the total electron/positron cascade y^ — 3^+ + 3^- can then be written as 

t 

ME, t, 9) J dt'e-^'^'''>^'-''^d^sye{E, t\ 9) - V,{E)y,{E, t, 9) + Cf{E, t, 9) 



(Al) 



j dt'e-r=(^)(*-'') j dE'j,,iE',E) 

E 



y,{E',t',9) + r,{E')y,{E',t',9) ~ cf{E',t', 



with an effective source term 



CfiE,t,9)= J dE'^^,{E',E)y^{E',t,9)+C:iE,t,9). 



(A2) 



(A3) 



a— e,7 



For tTe ^ 1 we can make the replacement Fg exp(— Fe(^ — t')) — > S{t — t') and Eq. (A2) reduces to 

djy,{E,t,9) 



Rl{E) 



dE'[r,{E')5{E^E')-j,e{E',E) 



y, {E\ t, 9) + r, [E')y,{E\ t, 9) - cf [e\ t 



(A4) 



We can further simplify Eq. ( A4 ) by introducing the mean inelasticity, 

.E' ^UE,E') 



= 1 - dE^ 



E r,{E) 



(A5) 



The inelasticity of ICS off CMB photons for electron/positron energies below about 100 TeV is small and we can 
hence approximate the differential interaction rate by ^ee{E' ,E) ~ Tf,{E')5{E — E'{1 — {x))). Using this in Eq. (A4) 
and taking the limit {x) ^ 1 we arrive at 



dly,{E,t,9) 
{x){E)Rl{E) 

Integrating this equation gives 



dE[ET,{E) 



dE 



„ djy, {E',t,9) 



ET^iE) J {x){E')Rl{E') 



y,{E, t, 9) + r,iE)y,iE, t, 9) - cf [e, t, 



ME,t,9) + r, [E)y, {E,t,9)- cf [E, t, 



(A6) 



(A7) 



We hence arrive at the diffusion term ^ with diffusion matrix ([7| for a regular magnetic field. 

If the coherence length As of the magnetic field is smaller than the distant to the source we can not neglect the 
spatial dependence of the diffusion velocity R£^. Generalizing to two angular variables 9 = {9i,92) in the plane 
orthogonal to the line-of-sight we start with 



3>± {E, t,9)=± -jrn^riL mey± (E, t, 9) ~ Ee iE)y± {E, t, 



dE' ( -j^e{E', E)y^{E', t, 9) + 7ee(£;', E)y±{E'., t, 



;C:iE,t, 



(A8) 
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where Ul is the direction of the Lorentz force projected into the 0i02-pla'ne. From here we arrive at 

t 



1 



RUE) 



dt'e-r=(^)(*-*')Tii(i)Ve \nL{t')Weye{E, t', 0)] - T,{E)y,{E, t, 6) + Cf{E, t, 0) 



Z 

J dt'e-r=(^)(*-*') J AE'^,,{E\E) 

E 



ye{E', t\ 0) + T,{E')y,{E\ t', 0) - CfiE', t\ 



(A9) 



For the evaluation of the second time integral in Eq. (A9) we can proceed as in the case of a regular magnetic 
field. However, in the first integral we have to account for the fluctuations of riL(t') over the inverse Compton 
scattering length. These will average to zero except for /Si < As and we hence substitute reexp(— re(i — t')) ^• 
min(l, ABre)(5(i — t'). Averaging over the orientation of the magnetic field can be accounted for by an additional 
factor 1/3. Proceeding now along the same steps as in the case of a regular field and replacing the angles 6*1/2 by 



spherical coordinates with radius we arrive at the diffusion term (10 1 with diffusion matrix (11) 



Appendix B: Diffusion-Cascade Equations 



We start from the Boltzmann equations (jlj and and define discrete values fJ"'' ~ AEiYe"'\Ei), Q^ i ~ 
AEiQe{Ei), etc. The combined effect of transitions and deflections within the cascade during a sufficiently small 
time-step At can be described by 



Y^{t + At)\ 
Yeit + At) 



Y^{t + At)\ 
Y,{t + At) 



-^WiAt) T^.iAt))^\YAt))^ 

\ (n) 



T^j(At) T,^{At) 

T^e{At) Tee (At) 



Y^{t) 

Yeit) 



3 

(n) 



Qe 


V 



Y^it) 

Yeit) 



(n-1) 



(n > 0) , 



The full cascade solution is then given by 



Y,{r 

Yeit'] 



in) 



m— j 



Y,{t) 

Yeit) 



(n— m) 



AtY^Bfit'-t) 



Qe 



(Bl) 
(B2) 



(B3) 



The 2n matrizes yt'™) and follow the recursive relation 

n 

(2PAt) = A^''' At) • (2^-1 At) , (B4) 

n 

S(")(2PAt) = S(")(2f-iAt) + ^y^('H2f-iAt) • S("-')(2P-iAt) , (B5) 

i=0 

where the non-zero initial conditions are A^^HAt) = r(At), A[f = diag(0, AtP.j) and B^°'>iAt) = 1. The matrices 

and S*^"^ are the familiar transfer matrices for electro-magnetic cascades in the presence of a source term. Using 
the recursion relations (B4) and (B5) we can efficiently calculate the matrices and via matrix-doubling |22j . 
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